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Abstract 

We show that in asymmetric heavy-ion collisions, especially off-central Cu+Au collisions, a siz- 
able strength of electric field directed from Au nucleus to Cu nucleus is generated in the overlapping 
region, because of the difference in the number of electric charges between the two nuclei. This 
electric field would induce an electric current in the matter created after the collision, which result 
in a dipole deformation of the charge distribution. The directed flow parameters vf 1 of charged 
particles turn out to be sensitive to the charge dipole and provide us with information about electric 
conductivity of the quark gluon plasma. 
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Introduction. — The quark gluon plasma (QGP), which consists of deconfined quarks 
and gluons, is expected to have filled the early universe We are now at the stage to 
study properties of the QGP experimentally through the relativistic heavy-ion collisions 
at Relativistic Heavy Ion Collider (RHIC) in BNL and Large Hadron Collider (LHC) in 
CERN. One of the most interesting observations is the very strong elliptic flow in off-central 
collisions, which indicates a very small ratio of shear viscosity to entropy density rj/s [25]. 
We are now trying to learn more detailed properties of the QGP by constraining transport 
coefficients such as shear viscosity, bulk viscosity, and charge diffusion constants. 

In this Letter, we propose a new way of estimating the electric conductivity a of the QGP 
via asymmetric nucleus-nucleus collisions at ultrarelativistic energies. Theoretically, lattice 
QCD simulations js-8] and perturbative QCD calculations {9] have been utilized to estimate 
electric conductivity of the QGP. So far, the estimated values of a have differed significantly 
from each other and experimental information is very much awaited. Very recently, asym- 
metric collisions between Copper (Cu) and Gold (Au) nuclei have been performed at RHIC 
and the PHENIX Collaboration reported their first results [loj]. We show that Cu+Au colli- 
sions can be useful for extracting the electric conductivity of the QGP. In off-central Cu+Au 
collisions, a substantial magnitude of electric field directed from a colliding Au nucleus to 
Cu nucleus is generated in the overlapping region. This happens only when colliding two 
nuclei carry the different number of electric charge. This electric field would induce a cur- 
rent in the matter created after the collision, resulting in a dipole deformation of the charge 
distribution in the medium. Later, the time evolution of the system is dominated by the 
strong radial flow, which is an outward collective motion of the medium. Henceforth the 
charge asymmetry formed in the early stage is frozen. Thus, we argue that charge-dependent 
directed flow of the observed hadrons are sensitive to the charge dipole formed at the early 
stage, which reflects the electric conductivity of the QGP. 

Conventionally the electric conductivity of the QGP is estimated from experiments via 
the Kubo formula The production rate of thermal dilepton is expressed by the electric 
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13] and its small frequency region is governed by the 



current-current correlation function 

transport peak jl^ . Thus one can estimate the electric conductivity through comparison of 
theoretical results with dilepton invariant mass spectra 15j . Compared to this method, the 
present approach is a rather direct one, in which the response of the matter to an applied 
electric field is directly quantified. 



2 



The effects of transient strong electromagnetic fields are under intensive discussions re 



cently, especially in the context of the chiral magnetic effect |16Nl8|. So far, there has 
been no experimental evidence that strong fields actually exist. Observation of a charge- 
dependent directed flow would also provide evidence that a strong electromagnetic field is 
actually created in heavy-ion collisions. 

Electric fields in Cu+Au collisions. — Here, we show that, in off-central collisions between 
copper and gold nuclei, a sizable strength of electric field is generated in the overlapping 
region of two nuclei. Because of the difference in the number of protons between the two 
nuclei, the generated electric field tends to the copper nucleus. The situation is different 



from the electromagnetic fields in the collisions of the same species of nuclei |19l . |20|. In 
symmetric collisions like Au+Au or Cu+Cu, the event-averaged electric field does not have 
a specific direction, although the magnitude of the electric fields generated in each event is 
considerably large {\eE\ ~ |eJ3| ~ 0{m 2 lv )). We have performed event-by-event calculations 
of the electromagnetic fields in Cu+Au collisions to show that there should be a significantly 
large copper-tended electric field. 

The electromagnetic fields are generated by the protons in nuclei. If we regard protons 
as point particles, the electric and magnetic fields at a position x and time t is written by 
the Lienard-Wiechert potentials, 



\e\E(t,x) = a EM ^ — 

\e\B(t,x) = a E M^ — 
» Ri 
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where R n = x — x n (t) with x n (t) the position vector of n-th proton, v n is the velocity vector 
of n-th proton, |e| is the electric charge of a proton, and «em is the fine structure constant. 
We have defined the origin of the spatial coordinate as the middle of the centers of the nuclei 
and x and y axes as in Fig. [TJ The summation is taken over all the protons in the colliding 
two nuclei. The positions of the protons inside a nucleus are sampled from the Woods-Saxon 



distribution with the standard parameters 



2l|. 



Figure [2] shows the event-averaged electric fields in Cu+Au collisions at impact parameter 
b = A fm. Each vector represents direction and magnitude of the electric field at that point. 
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FIG. 1: (Color online) Transverse plane of off-central Cu+Au collisions with impact parameter b. 
The left (right) circle indicates the edge of the Cu (Au) nucleus. 

We find that the electric field in the central region of the overlapping area has a specific 
tendency to go from Au to Cu. Although the direction of electric fields fluctuates on an 
event-by-event basis because of the fluctuation in the proton positions inside colliding nuclei, 
the direction is correlated with the reaction plane for asymmetric collisions. The magnitude 
of the electric fields is as large (\eE\ ~ 0{m? 1T )) as the electric and magnetic fields in Au+Au 
collisions at the same collision energy. 

We have also calculated the time dependence of the averaged electric fields as shown 
in Fig. El The strength of the fields decays as the spectators fly away. Nevertheless, it is 
notable that even at t = 1 fm/c the electric field is considerably larger than the so-called 
"critical field" for electrons, \e\B c = \e\E c = m\. 

Electric dipole of the plasma and charge- dependent directed flow. — The strong electric 
field toward the Cu nucleus at the early stage would induce an electric current in the medium 
that consists of the QGP after the thermalization time. As a result, the charge distribution 
would be modified and a charge dipole would be formed. One can expect that the dipole-like 
deformation of the charge distribution at the early stage would also be present in the observed 
charge distribution. This is because the electromagnetic charge is an exactly conserved 
quantity and an inhomogeneity relaxation of a conserved charge density takes a long time. 
Once a radial flow starts, the medium expands rapidly and the charge dipole created at the 
early stage would be frozen. Thus, we can reasonably assume that the dipole-deformation 
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FIG. 2: (Color online) Event-averaged electric field in the transverse plane in off-central Cu+Au 
collisions at t = (the collision time) with impact parameter b = 4 fm at ^/snn = 200 GeV. 
Vectors are shown only in \y\ < 6 fm. The average is taken over 10 4 events. 
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FIG. 3: (Color online) Event average of the time evolution of electric fields in off-central Cu+Au 
collisions (6 = 4 fm) at x = 0. The average is taken over 10 4 events. 

in the plasma remains intact in the observed charge distribution. 

The azimuthal angle distribution of the net charge, to the leading order in the multipole 
expansion, can be written as 

d (N+ d ~ — (0) = (JV+ - N-) (1 + 24 cos 0) , (3) 
where the azimuthal angle <j) is measured from the x-axis and N± is defined as the angle 



average of the number distribution, 



N ± EE / (4) 



The dipole deformation of the medium is quantified by the value of d e . We assume that the 
azimuthal distribution of the total number of particles is still written by the vx without the 
effect of the electromagnetic fields as 

d{N+ ^ N ~ ] = (#+ + N-) (1 + 2v, cos0) , (5) 

since electromagnetic fields are not expected to change the bulk flow significantly. From 
Eqs. ([3]) and ([5]), the distribution of charged particles can be written as 



1 H ^ 2 (vi + Ad e ) cos i 



(6) 

= N + {1 + 2 [vx + A(d e - Vl )\ cos0 + O [(Ad e ) 2 ] } , 

where we have defined the charge asymmetry parameter A = (N + — N_)/(N + + N_). 
Similarly, 

^- = N_{l + 2[ Vl - A(d e - Vl )} cos0 + O [(Ad e ) 2 } } . (7) 

Thus, the directed-flow coefficients V\ for positively and negatively charged particles are 
written as 

vf = vx± Ad' e , (8) 

where we have defined d' e = d e — Vx- The values vf are linear functions of A and their slopes 
are given by the dipole-like deformation parameter d' e , which is written as 

4 = ^—^ J rdrd( P cos jl{r,4>) -3 ei E=B=o( r > < t > ) ' ( 9 ) 

where je(r,4>)(j®g_g_ Q (r,(f))) is the transverse charge density in the presence (absence) of 
electromagnetic fields. 

Estimate of the charge- dependent directed flow. — Let us make an order-of-magnitude 
estimate of the value of the charge-dependent directed flow parameter Ad' e . For that purpose, 
we first roughly evaluate the total charge that is transfered from the gold-side to copper-side 
in the presence of an electric field. The total charge Q transfered across a plane S from t — 
to r is written as 

Q= f dt [ J-dS= [ dt [ oE-dS, (10) 
Jo J s Jo Js 
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where we have used the constitutive relation J = aE with a the electric conductivity. Let 
S be the plane which includes the origin and is perpendicular to the line which connects the 
centers of the two colliding nuclei at t — 0, the moment two nuclei contact. Neglecting the 
space-time dependence of a, Q is rewritten as 

Q~ctt [ E-dS. (11) 
Js 

The integral in Eq. (TTTj) is just the total electric flux that goes through the plane S. Hence, 
the total transfered charge Q is roughly given by 

E-dS- Za "" Zcu M, (12) 
s 2 e 

where Z^ n and Zq u are the number of protons in two nuclei, and e is the dielectric constant 
of the QGP. 

According to lattice QCD simulations, the electric conductivity of the QGP is estimated 



as 



20] 



a ~ B CemT, C E M = ^ej, (13) 

/ 

where the sum in the electromagnetic vertex factor is taken over the flavors and B is a 
coefficient that depend on the calculations. If we consider u, d, and s quarks, Cem — 
87to;em/3. The value of the coefficient B differs among calculations: B ~ 0.4 in Refs. [7J, |8| 
and B ~ 7 in Ref. jfj]. On the other hand, perturbative QCD calculations predicts a ~ QT/e 2 
(9], which is much larger than the values from lattice QCD simulations. 

As for r, we take the time scale that the radial flow starts, r ~ 1 fm/c. If we take typical 
values for the other parameters, T ~ 200 MeV, e ~ 1, and B ~ 1, the total transfered 
charge is estimated as 

Q ~ iiC EM i t — 

~ 1 • -«em ■ 200 MeV- 1 fm/c- 25|e| ( 14 ) 
3 

~ 1.7 |e|. 

Now we can roughly evaluate Ad' e . Let us choose the events in which the numbers of 
positive and negative hadrons are equal, N + = N-, and assume that n charges have been 
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transferee! by the electric field. Then, the number n can be written as 22] 

1 r /2 d(iV+ - AL) 



n = "2 



tt/2 a<P (15) 

= -2A< (A\ + At) . 

Therefore, the directed flow parameter Ad' e is written by n as 

Ad' e = (16) 

where N tot = 2tt (iV + + iV_) is the total number of charged particles. The number n is 
related to the total transfered charge roughly as n ~ Q/|e|. Therefore, Ad' e , the charge- 
dependent part of the directed flow parameter, and the electric conductivity of the plasma 
are parametrically related as 

Ad> e ~-^p- I E-dS. (17) 

- ,v totp| Js 

If one takes Atot ~ 10 3 and n ~ 1 (Eq. fTH|) ). the order of magnitude of the directed- flow 
parameter is estimated as 

Ad' e ~ -10~ 3 . (18) 

This value would be within experimental reach. Note that the value (1181) is negative since 
the electric field tends the Cu nucleus. 

Let us comment on potential uncertainties in the estimate above. It would be possible 
that the charge dipole formed at the early stage can be obscured in the later stages, namely 
the hydrodynamic evolution and hadronic collisions. In order to quantify these effects, we 
have to calculate the time evolution of the charge density under an electric field. Hence, 
it is necessary to use a hydrodynamic model which includes electromagnetic fields and a 
hadronic afterburner. One should also consider the effects of fluctuations of the generated 
electric fields on an event-by-event basis, although the fields have a tendency to direct from 
Au to Cu nucleus on average. The charge dipole could be smeared out by the fluctuations. 
Event-by-event simulations are necessary to consider the effect of such fluctuations. Finally, 
although we have assumed the dielectric constant is a constant, it can in general depend on 
frequency and wave length. Consideration of such effect is left as a future work. 

Summary. — We have pointed out that, in Cu+Au collisions, a sizable strength of electric 
field directed from Au to Cu nucleus is generated in the overlapping region. We have shown 
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this by performing event-by-event numerical calculation of the produced electromagnetic 
fields. We have also pointed out that the electric field would induce an electric current in 
the matter created after the collision and it would result in a dipole deformation of the 
charge distribution in the medium. We have shown that the charge- dependent directed flow 
of hadrons is sensitive to the charge dipole in the medium and is useful in estimating the 
electric conductivity of the QGP. 
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